Influence of long-range correlated quenched disorder on the 
adsorption of long flexible polymer chains on a wall. 

Z. Usatenko 1,2 , A.Ciach 2 
Institute for Condensed Matter Physics, 
National Academy of Sciences of Ukraine, 
1 Svientsitskii Str., UA-79011 Lviv, Ukraine * and 
2 Institute of Physical Chemistry of the Polish Academy of Sciences, Warsaw, 01-224, Poland 

Abstract 

The process of adsorption on a planar wall of long-flexible polymer chains in the medium with 
quenched long-range correlated disorder is investigated. We focus on the case of correlations 
between defects or impurities that decay according to the power-low x~ a for large distances x, 
where x = (r, z). Field theoretical approach in d = 4 — e and directly in d = 3 dimensions up 
to one- loop order for the semi- infinite |t/>| 4 m- vector model (in the limit m — > 0) with a planar 
boundary is used. The whole set of surface critical exponents at the adsorption threshold T = T a , 
which separates the nonadsorbed region from the adsorbed one is obtained. Moreover, we calculate 
the crossover critical exponent <3? and the set of exponents associated with them. We perform 
calculations in a double e = 4 — d and 5 = 4 — a expansion and also for a fixed dimension d = 3, 
up to one-loop order for different values of the correlation parameter 2 < a < 3. 

The obtained results indicate that for the systems with long-range correlated quenched disorder 
the new set of surface critical exponents arises. All the surface critical exponents depend on 
a. Hence, the presence of long-range correlated disorder influences the process of adsorption of 
long-flexible polymer chains on a wall in a significant way. 
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I. INTRODUCTION 



Universal properties of long-flexible polymer chains change when small amount of long- 
range correlated quenched disorder is introduced into an infinite medium [1, 2]. The corre- 
lated defects (i.e. regions that cannot be occupied by the chain) may occur, for example, 
in a porous medium or in a disordered sponge-like structure formed by lipid membranes 
in biological systems. Intuitively, if the correlations between the defects and/or impurities 
decay sufficiently slowly, then the chain has to go around large correlated regions, and ef- 
fectively occupies larger space, with the defects contained inside the region occupied by the 
coil. As a result, the polymer swells. If, however, the range of correlations is very large, 
then the polymer may be trapped between the walls of defects (i.e. the probability of going 
beyond the defected region is low), and this may lead to a collapsed state. These heuristic 
arguments suggest that the polymer either swells or collapses, depending on the range of 
correlations between the defects or impurities. Indeed, recent results agree with the intu- 
itive expectations [1, 2]. For different ranges of correlations (i.e. different values of a for the 
power-law decay of correlations l/x a ) the swelling of the polymer is described by different 
dependencies of the radius of gyration on the number of monomers. Finally, for a < 2.3 
[1, 2] a first-order transition to a collapsed state was found [1, 2}. 

Motivated by the above results we focus our attention on the effect of the presence of 
a small amount of long-range correlated quenched disorder in the bulk on the adsorption 
of long-flexible polymer chains on a planar surface forming the system boundary. In real 
systems different kinds of defects and impurities may be localized inside the bulk or at 
the boundary. As was found in Ref. [15], introducing into the system short-range correlated 
random quenched surface disorder is irrelevant for critical behavior, but long-range correlated 
quenched surface disorder with g(r) ~ l/r a can be relevant if a < d — 1, and is irrelevant 
if a > d — 1. The question how the adsorption phenomena of long flexible polymer chains 
depend on the presence of long-range-correlated quenched disorder in the bulk remains open, 
however. New universality class characterizing the polymer in the presence of long-range 
correlated disorder indicates that the critical exponents describing the properties of the 
polymer chain near the wall should assume different values than in the pure system. The 
purpose of this work is a determination of the surface critical exponents to first order in 
perturbation expansion, in order to gain information about a qualitative change of adsorption 
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of the chains when the range of correlations between the defects in the bulk increases (i.e., 
a is decreased). 

Long flexible polymer chains in a good solvent are perfectly described by a model of 
self-avoiding walks (SAW) on a regular lattice [16]. Their scaling properties in the limit 
of an infinite number of steps N may be derived by a formal m — > limit of the 0(m) 
vector model at its critical point [3] . The average square end-to-end distance, the number of 
configurations with one end fixed and with both ends fixed at the distance x = \J (xa — xb) 2 
exhibit the following asymptotic behavior in the limit N — > oo 

<R 2 >~N 2u , Z N ~ q N W\ Z N {x) ~q N N-( 2 ~ a \ (1) 

respectively, v, 7 and a are the universal correlation length, susceptibility and specific 
heat critical exponents for the m = model, d is the space dimensionality, q is a non 
universal fugacity. 1/N plays a role of a critical parameter analogous to the reduced critical 
temperature in magnetic systems. 

When the polymer solution is in contact with a solid substrate (or with vapor), then the 
monomers interact with the surface (or their chemical potential at the interface is different 
than in the bulk). At sufficiently low temperatures, T < T a , the attraction between the 
monomers and the surface leads to an adsorbed state, where a finite fraction of the monomers 
is attached to the system boundary. Deviation from the adsorption threshold, c oc (T — 
T a )/T a ), changes sign at the transition between the adsorbed (c < 0) and the nonadsorbed 
state (c > 0) and it plays a role of a second critical parameter. The adsorption threshold for 
infinite chains, where 1/N — > and c — > is a multicritical phenomenon. We shall assume 
that the solution of polymer chains is sufficiently dilute, so that interchain interactions and 
overlapping between different chains can be neglected, and it is sufficient to consider surface 
effects for configurations of a single chain. For pure solvents the investigation of adsorption 
phenomena of long- flexible polymer chains on a surface was a subject of a series of works 
(for the sake of brevity we notice only few of them [3, 4, 5, 6, 7, 8, 9, 10, 11]). The polymer 
adsorption on a wall in the limit of an infinite chain is closely related to surface critical 
phenomena in the m - vector model of a magnet in a semi-infinite geometry in the limit 
m — > [3, 17, 18]. Based on the above analogy, Eisenriegler and co-workers [4, 5, 6] described 
the scaling properties of long chains near a wall on the basis of the results of the field 
theory developed for semi-infinite magnetic systems in Ref.[12, 13, 14]. Surface multicritical 
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phenomena in dilute polymer systems (at T = T a and iV — > oo) correspond to the special 
transition in semiinfinte magnets. The special transition (c = 0) is characterized by one 
additional independent surface critical exponent r)\\, which characterizes critical correlations 
in directions parallel to the surface. The whole set of the other surface critical exponents 
can be obtained on the basis of r]\\ and the bulk critical exponents v and r\ with the help 
of surface scaling relations. The crossover critical exponent $ characterizes the crossover 
behavior between the special and ordinary transitions (c^O). The latter exponent is related 
to the length scale [5, 6], 

tc~\c\- v/ ; (2) 

associated with the parameter c. In the polymer problem the length £ c can be interpreted 
as the distance from the surface up to which the properties of the polymer depend on the 
value of c, not only on its sign. The remaining, bulk length scales are the average end-to-end 
distance £ R = V< R 2 > ~ N v and the microscopic length / - the effective monomer linear 
dimension. Near the multicritical point the only relevant lengths are £r — > oo and £ c — > oo, 
and the properties of the system depend on the ratio £r/£ c - In the asymptotic scaling regime 
the universal physical quantities X(N, c) and Y(z; N, c) assume the scaling forms 

X(N, c) = N a -X^ R /^ C ), Y(z; N, c) = N a -Yl(z/^ (3) 

where X± and Y± denote the scaling functions with the subscripts + and — corresponding 
to c > and c < respectively. The characteristic length ratio is (Cr/^c)*^ ~ M^*) where 
cN® is the standard scaling variable [4] . The exponents ax and ay assume different values 
for different quantities X and Y. Let us first consider the mean square end-to-end distance 
for one end attached to the surface and the other one free. In the semi-infinite system the 
translational invariance is broken, and the parallel < R?, > and perpendicular < R\ > parts 
of the average end-to-end distance < R 2 >=< R\ + R 2 > should be distinguished. For 
< R\ > l l 2 the exponent in the scaling form (3) is ax = v and the corresponding scaling 
functions assume the form ~ const for c > and ~ 1/y for c < 0, where y = £j?/£ c [4]. 
Thus, for the adsorbed state and for iV — > oo the length associated with c describes the 
thickness £ of the adsorbed layer, 

£ =< Rl >V*„ £ c c < 0. (4) 
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This thickness diverges for c = 1/JV = and for finite negative values of c remains finite for 
an infinite chain. For c > the asymptotic behavior of the mean distance of the free end 
from the other end attached to the surface is 

< R\ > 1/2 ~ N v c > 0, (5) 

i.e it has the same asymptotic behavior as in the bulk. The asymptotic scaling form of 
< Rl > 1/2 for c < is < Rj > 1/2 ~ | c |( I ' d-1 - I ')/*jV I ' d_1 ) where v d ~ l is the correlation exponent 
in d — 1 dimensions. For c > the scaling form of < i? 2 > x l 2 is given by Eq. (5), i.e. it is 
also the same as in the bulk. 

For the fraction of monomers at the surface, Ni/N, the following asymptotic behavior 
has been found for iV — > oo [4, 6] , 

' | c |(i-*)/« if c<0 

N^N ~ j JV*" 1 if c = • (6) 
_ (cN)- 1 if c> 

Hence, for N — > oo and for finite, negative values of c, N±/N is finite, but for c > 
Nx/N — > for A" — > oo. The thickness of the adsorbed layer is closely related to the fraction 
of monomers at the surface Ni/N [4, 6], since the more monomers are fixed at the wall, the 
smaller the region occupied by the remaining monomers. In particular, for weakly adsorbed 
phase (c < and |c| < 1) we find NJN ~ f-t 1 "*)/". 

The scaling behavior is also obeyed by the mean number of the free ends in the layer 
between z and z + dz, which is proportional to the partition function of a chain with one 
end fixed at = (r^, z) and the other end free, Z N (z). The density of monomers in a layer 
at the distance z from the wall to which one end of the polymer is attached, Mj^(z) scales 
according to Eq. (3) as well. For the above quantities the exponent ay in (3) is 7 — 1 and 
71 — v respectively. Short-distance behavior (/ <^ z <C £_r) of the two quantities right at the 
threshold (c = 0) is 

Z N (z) ~ Z ln-*)I"N^- 1 (7) 

and 

Mfa) ~ z -i+(i-*)/-AT7i-i+*. ( 8 ) 

71 as well as the whole set of the surface critical exponents can be obtained from r)\\ and $ 
through scaling relations (see the Appendix). The remaining quantities characterizing the 
adsorption process are described in detail in Ref.[4, 6]. 



Taking into account the results of ref. [1, 2] we conclude that for c > the polymer with 
one end attached to the surface swells as in the bulk when a decreases (see (5)). However, 
in order to determine the effect of the long-range correlated disorder on the adsorption of 
the polymer right at the threshold (see (7) and (8)) or in the crossover region (see (2), (4) 
and (6)) it is necessary to find the dependence of the surface critical exponents on a. 

In the next section the model is briefly described. In sec. Ill the surface multicritical 
behavior of the system with long-range correlated disorder is outlined. The results of sec. Ill 
enable us to obtain in sec. IV the surface critical exponents to first order in the perturbation 
expansion. Final section contains a brief discussion of the results. 

II. THE MODEL 

When a disorder is introduced into an infinite magnetic system, the Landau-Ginzburg- 
Wilson Hamiltonian [3] assumes the form 

H = j v d d x[\\ V0| 2 +^(^ + 5r(x))|0| 2 +^o(^) 2 ], (9) 

where <f>(x) is an m-vector field with the components <f>i(x), i = l,...,m. Here /z 2 , is 
the "bare mass", which in the case of a magnet corresponds to the reduced temperature. 
The inhomogeneities in the system cause local deviations from the average value of the 
transition temperature, and 5r(x) represents the quenched random-temperature disorder, 
with < 5r(x) >= and 

i<5r(x)Mx)>=s(|x|), (10 ) 

where angular brackets < ... > denote configurational averaging over quenched disorder. 
Following Refs.fl, 2] we assume that the pair correlation function g(\ x |) falls off with the 
distance as 

(11) 

for large x = (r, z), where a is a constant and x —\ x |. The Fourier-transform g(k) of g(x) 
for small k is 

~g{k) ~ u + w | k \ a ' d . (12) 

This corresponds to the so-called long-range-correlated "random-temperature" disorder. In 
the case of random uncorrelated point-like (or short-range-correlated) disorder the site- 
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occupation correlation function is g(x) ~ S(x) and its Fourier-transform assumes the simple 
form 

g(k) ~ u Q . (13) 

Applying the replica method in order to average the free energy over different configu- 
rations of the quenched disorder, it is possible to construct an effective Hamiltonian of the 
|0| 4 m -vector model with a long-range-correlated disorder 

n r 1 1.1. 
Heff = E//^ I V0 a | 2 + ^v (& 

n . 

- E / (fxKfx'^d X! - x[ DfcWfH*!)- (14) 

a,(9=l 

Here Greek indices denote replicas, and the replica limit n — » is implied. If a > d, then 
the wq term is irrelevant. This corresponds to random uncorrelated point-like disorder (or 
short-range-correlated random disorder). As noticed by Kim [20], in this case in the limit 
m, n — > both vo and uq terms are of the same symmetry. It indicates that a weak quenched 
uncorrelated disorder is irrelevant for SAWs [21]. If, on the other hand, a < d, the term 
Wok a ~ d is relevant for the critical behavior at k — > 0, and the long-range-correlated disorder 
is relevant for SAWs (see [1, 2]). The limit m — > of this model can be interpreted as a 
model of long-flexible polymer chains in a disordered medium [1, 2]. 

The presence of a hard wall leads to a modification of the interactions in the near-surface 
layer. Thus, in the semi-infinite system there should be additional, surface contribution to 
the Hamiltonian. The effective Hamiltonian of the semi- infinite |0| 4 m- vector model with 
long-range-correlated disorder in this case is 



n r 1 1^1^ 
Heff = E/ y ^ 1 V0a |2 +2^ + 4T o(0 ' )2 

- E / d^d^gd Xl - Dg^)^^) + ^ E / ^-V0l(r, z = 0),(15) 

a,p=l a=l 

where Co describes the surface-enhancement of interactions. In the polymer analog Co oc 
(T — T a )/T a , as already noted in the introduction. The surface introduces an anisotropy 
into the problem, and directions parallel and perpendicular to the surface are no longer 
equivalent. In accordance with the fact that we have to deal with semi-infinite geometry 
(x = (r, z > 0)), only parallel Fourier transforms in d — 1 dimensions will be performed. 

7 



The parallel Fourier transform g(q,z) of (10) is 

g(q,z)=w d q^z^^Ka^+iiqz), (16) 

where K a-d+i (qz) is the modified Bessel function and q —\ q |, where q is a d — 1 dimensional 
vector. In the case of small q and z we obtain the relation 

g(q, z)~u + w q a - d+1 + w' z d - a -\ (17) 

which agrees with the predictions obtained in [15]. We concentrate our attention on the case 
a < 3 for d = 3, for which the long-range correlated disorder in the bulk is relevant. In the 
general case of arbitrary z (from z = on the wall to z — > oo) we must take into account 
the Fourier transform g(q, z) of the form (16). 



III. SURFACE CRITICAL BEHAVIOR NEAR THE MULTICRITICAL POINT 



c = c„ 



A. Normalization conditions 

The correlation function which involves N fields 0(x;) at distinct points Xj(l < % < N) 
in the bulk, M fields <j>(rj,z — 0) = (j> s ( r j) at distinct points on the wall with parallel 
coordinates rj(l < j < M), and L\ insertions of the surface operator |^(Rj) at points R; 
with 1 < I < Li, has the form 

{r,>, {r}) =< n ^) n n ^'( R ') >. 



i=l j=l z=i " 

where < ... > denotes averaging with the Boltzmann factor, in which the Hamiltonian is 
given in Eq.(15). The corresponding full free propagator in the mixed pz representation is 
given by [14] 



e -K0\z-z'\ _ °0 ~ K c -K n (z+Z') 



(19) 



C + Kq 

where n = \Jp 2 + /Xq with p being the value of the parallel momentum p associated with 
the d — 1 translationally invariant directions in the system. 

There are two special cases: a) when two ends of the polymer are attached to the wall (in 
such a case we have to deal with a calculation of two point correlation function G^ ^ (r, z = 



0; r', z' = 0)), and b) when one end of the polymer is unrestricted in the bulk and the other 
one is attached to the wall (G^^^x; r' , z' = 0)). In order to obtain the universal surface 
critical exponents, characterizing the adsorption on the wall of long-flexible polymer chains 
inserted into the medium with long-range-correlated quenched disorder, it is sufficient to 
consider the correlation function of two surface fields G^ ,2 \r, z = 0;r',z' = 0) (see [19]). 
The universal surface critical exponents for such systems depend on the dimensionality of 
space d, the number of order parameter components m{m — > 0) and the range of the disorder 
correlations, i.e. on a. 

In the theory of semi-infinite systems the bulk field 0(x) and the surface field s (r) should 
be reparameterized by different uv-finite renormalization factors Z^{u,v,w) and Zi(u,v,w) 
[12, 19], 

0(e) = Z^ 2 <f> R {x) and (f> s (r) = Z l J 2 Z\ /2 (f) SyR (r). 

Introducing the additional surface operator insertions |0 2 (R;) requires additional specific 
renormalization factor Z^ s 

The renormalized correlation function involving N bulk, M surface fields and L x surface 
operators |0^(R;) can be written as 

4 Wl W, W,c) = ^^/ 2 Zr M/2 ^G^)(p ;/ , , Mo ,t;o, W o,Co). (20) 

It should be mentioned that the typical bulk short-distance singularities, which are 
present in the correlation function G^ ,2 \ can be subtracted after performing the mass renor- 
malization. For distinguished parallel and perpendicular directions we obtain: 

^ = /i 2 -4°W) + 4 0) / 2 (A (21) 

where 



t^- l (v v \ .(Q) _ (22) 



'd-a\ ' 

and 



^ 2 ' = (2^/^7 ^ 
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According to the above mentioned notations, we have only two coupling constants, V = 
t>o — Mo and wo in the effective Hamiltonian (we keep notation vq for Vq). 

The renormalized coupling constants v, w are fixed via the standard normalization con- 
ditions of the infinite-volume theory [2] 

^- d v = r ( ^ v {{q};fi 2 ,v,w)\ q= o, 

^ = ri 4 U{g};/i 2 ,^)| 9= o, ( 24 ) 

where T^ R ^ and T^ R w are the v and w term symmetry contributions to the four-point vertex 
function. To the present accuracy of calculation at one-loop order, the vertex renormalization 
gives: v = Vofi d ~ 4 and w = Wofi a ~~ 4: . 

In order to remove the short-distance singularities of the correlation function G^ 0,2 -*, lo- 
cated in the vicinity of the surface, the surface-enhancement shift Co = c + 5c is required. In 
accordance with this, a new normalization condition should be introduced for the surface- 
enhancement shift 5c and the surface renormalization factor Z\. By analogy with magnetic 
systems [19, 24, 25], the renormalized surface two-point correlation function in our case is 
normalized in such a manner [19] that at zero external momentum it should coincide with 
the lowest order perturbation expansion of the surface susceptibility X||(p) = G^ ' 2 ' (p) 

G (0 ' 2) Go; Ho,v ,w ,c ) = 1 +O(v ,w ). (25) 

c + \JP 2 + n 

Thus, we obtain the necessary surface normalization condition, 

G% 2 \0;^v,w,c) = ^-, (26) 
{1 ~\~ c 

and for the first derivative with respect to p 2 we have 

dG { x' 2 \p;Li,v,w,c) 



dp 2 



2/x(/x + cf 

p=0 



(27) 



Eq. (26) defines the required surface-enhancement shift 5c and shows that the surface sus- 
ceptibility diverges at /i — c — 0. This point corresponds to the multicritical point (//q c , eg), 
at which the adsorption threshold takes place (it corresponds to the special transition). 
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From the normalization condition of Eq. (27) and the expression for the renormalized 
correlation function of Eq. (20), we can find the renormalization factor Z\\ = Z\Z^ from the 
relation 



p 2 =0 



p^o p dp 



(28) 



The normalization condition for the correlation function G^' 2,l \ with the insertion of the 

1 



surface operator \tp 2 s , 



/-»(0.2,1)/ \ 

Gr J {p;fi,u,v,c) 



P=0 (/i + c) 2 

gives a possibility of obtaining the renormalization factor Zxi from 



[Z 4 



dc 



(29) 



(30) 



co=co(c,fj,,u,v) 

Eq.(29) follows from the fact that the bare correlation function G^°' 2,1 \0; // , u 0: v 0: c ) may- 
be written as a derivative —■^G^°' 2 \0; Hq,u ,v , c ). 



B. The Callan-Symanzik equations 



Asymptotically close to the critical point the renormalized correlation functions G R 
satisfy the corresponding homogeneous Callan-Symanzik (CS) equations [19, 24] 

d d d N + M 



(N,M) 



M 



G ( R N ' M \0; f i,v,w,c) = 0, 



where the /^-functions are (3 v (v,w) = //^ v, (3 w (v,w) = piJ^ 
the exponents r\ and rfy are 



LR 



(31) 



sp 



9 1 7 I 



(32) 



and where LR is the long-range fixed point. It should be mentioned that up to one loop order 
in the e- and 5- expansion the LR is located in the region of irrelevant disorder a > 3, and 
up to two loop order the LR stable fixed point is found after performing the Borel-Chisholm 
resummation [1, 2]. 
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The simple scaling dimensional analysis of and of the mass dependence of the Z 

factors, allows to express the surface correlation exponent r/jj p as 

V\\ P = VT + V- (33) 

From Eqs. (28), (32) and (33), we obtain for the surface correlation exponent rj^ p the following 
expression 



LR 



n . . <91nZ||(t>, w) n . . d\nZ\\(v, w) 



(34) 



LR 



ov ow 
where the (3 functions are [2] 

Pv(v,w) = -v + v 2 - (3/i(a) - f 2 {a))vw - 

1)UJ 

P*(v,w) = -{A - a)w - (h{a) - f 2 {a))w 2 + — + ... (35) 

In the above equation the renormalized coupling constants v and w are normalized in a 
standard fashion, so that 

v = -vh, w = -wh, 

and the integral I\ in the case of d = 3 is equal to l/8n and in the case of d = 4 — e it is 
Ii = 2~ d ir~ d / 2 T(e/2). The coefficients fi(a), expressed via the one-loop integrals [2, 26] are 
given by 

(a-2)(a-4) 



2sin(7ra/2) ' 



(a-2)(a-3)(a-4) 
/2(a) = 4 87 rsim>(a/2-l)) " (36) 

C. crossover between the adsorbed and nonadsorbed states 

As already discussed in the introduction, it is particularly interesting to investigate the 
adsorption threshold and the crossover behavior between the adsorbed and the nonadsorbed 
states, where the distribution of monomers in the near-surface region changes character. 
In order to investigate the crossover behavior from the nonadsorbed region, c > eg, to the 
adsorbed one, c < eg, let us consider a small deviation Ac = c — Cq from the multicritical 
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point. The power series expansion of the bare correlation functions G^ N ' M \p; /i , v , w , c ) 
in terms of this small deviation Aco has the form 

G^ M \p;^v ,w ,c ) = £ { -^^G^ M ^(p^ ,v ,w ,c a ). (37) 

L 1= L ^ 

Taking into account Eq.(20), we can rewrite the right-hand side of Eq.(37) in terms of the 
renormalized correlation functions and renormalized variable Ac = [Z^^v, w)]~ 1 Acq. In this 
way we obtain 

Z -<» + »)/2 (Zi) _M/ 2G ( N ,M, (p . ^ Wo Co) 
L 1= L l ! 

The above equation determines in a straightforward fashion the corresponding renormalized 
correlation functions in the vicinity of the multicritical point (/iQ C , eg), 

G%> M \p^,v,w, Ac) = ^-^^(ZO-^G^Cp^o^o^o^o). (39) 

These correlation functions depend on the dimensionless variable c = Ac//i. The correlation 
functions G^' M ^(p; v, w, Ac) satisfy the CS equations (31)(see also Ref. [19, 24]) with the 
additional surface related term —[1 + r] £ (v , w)]c-^=, where 



( \ 9 
rjc{v,w) = fi— 



(40) 

LR 



should be calculated at the LR stable fixed point. 

The asymptotic scaling critical behavior of the correlation functions can be obtained 
through a detailed analysis of the CS equations, as was proposed in Ref. [27, 28] and 
employed in the case of the semi-infinite systems in [19, 29]. Taking into account the 
scaling form of the renormalization factor of Eq. (30) and the relation /i ~ t u , where 
t = (T — T c )/T c is the reduced critical temperature in magnetic systems, we obtain for Ac 
and for the scaling variable c the following asymptotic forms 

Ac ~ aT^^Aco, Ac ~ r-^^'^'Aco (41) 

and 

c-^+^'^Aco, c~r-*Ac , (42) 

where 

$ = v (l + n 5 {v*,w*)) (43) 
13 



is the surface crossover critical exponent. Eq. (42) explains the physical meaning of the 
surface crossover exponent as a value which characterizes the measure of deviation from the 
multicritical point. 

Taking into account the above mentioned results, we obtain from the CS equation the 
following asymptotic scaling form of the surface correlation function G^°' 2 \ 



point. The knowledge of r\c gives access to the calculation of the critical exponents ol\ and 
ol\\ of the layer and local specific heats via the usual scaling relations [14] 

a 1 = a + is-l + $ = l-is(d-2- rj 5 ), a\\ = a + v - 2 + 2$ = -v[d - 3 - 2rf B ]. (45) 

IV. THE PERTURBATION EXPANSION FOR THE SURFACE CRITICAL EX- 
PONENTS 

Applying the field-theoretical renormalization group (RG) approach we perform calcula- 
tions in a double expansion in e = 4 — d and in 5 = 4 — a up to the linear approximation, as 
was proposed by Weinrib and Halperin [30] for infinite systems. Thus, after performing the 
integration of the corresponding Feynman integrals in the renormalized two point correlation 
function G^ 0,2 \ we obtain at the first order of the perturbation theory the following result 
for the renormalization factor Z\\ 



G^ 2 \p; /i , v , w , c ) ~ //-^g£' 2) (-; 1, v*, w*, ^Ac ) 



~t-^G( P t-»;1,t-*Ac ), 



(44) 



where = u(l — r}\\), and r/y 



rj S i + r] are the surface exponents at the multicritical 




(46) 



where 



e l + o 



2F(S - e) 



x 




(47) 
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Combining the renormalization factor Z\\ together with the corresponding /^-functions de- 
rived in Ref. [1, 2], we obtain for the surface critical exponent r/| p the result 

sv v e ™5 . 

"ll'--4(rT7) + 2cos(^-.)) g( ^ ) - (48) 
The above mentioned surface critical exponent vf^ in the case of e, <5 - expansion can be 
calculated formally at the corresponding fixed point v* = ^rrjy, w* = — ^zfy^ obtained in 
the first order of e, 5 - expansion in [1]. 

In the special case of three spatial dimensions d = 3 and for arbitrary a the renormaliza- 
tion factor Z\\ at the one-loop order is given by 

Z|| - 1 + 8 + 8sin(vra/2)' (9) 
where we have introduced the function (7(a) by 

a-3 / 1 V 2 / , /o V 2 M 



.MfEffl + Iffl) (50) 

From the Eq.(50) it is easy to see that in the case of short-range correlated (or uncorrelated) 
disorder, i.e. for a = d = 3, the above mentioned function reduces to g(a) = 1, and both 
v and w terms are of the same symmetry. It confirms that a short-range correlated (or 
uncorrelated) disorder is irrelevant for SAWs. 

Combining the renormalization factor Zn together with the one-loop pieces of the f3- 
functions, according to Eq.(34), we finally obtain the following expression for the surface 
critical exponent r]* p , 

_ _v w (4 - a)g(a) 
^ ~ 8 8 sin(7ra/2) ' (5ij 
Similarly, for the renormalization factor Z^ we obtain at the one-loop order 

where h(a) is a combination of the Appell hypergeometric functions of two variables 

F 1 [a,b 1 ,b 2 ,c,x,y], 

h(a) = 2"- 4 F 1 [l, 1,^,2, -1,-1/2] + 

+ ^^[^ 1, ^ -1, -2] - t^S^t. (53) 

5 — a 2 2 2 (7 — a) (6 — a) 
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Finally, for the exponent f] 5 we obtain 

v /r n lx w{& — a) ., . , q(a). , . 

' te = -2< ta2 -i»-2^(wi)< fc ' a »-¥»- < 54 > 

In the case of the short-range correlated (or uncorrelated) disorder for the function h(a) at 
a = d = 3 we obtain ln2. 

The above values of the surface critical exponents 7/jj p and r] 5 should be calculated at the 
long-range (LR) stable fixed point obtained in Ref. [2, 26] for different fixed values of the 
correlation parameter, 2 < a < 3. The other surface critical exponents can be calculated on 
the basis of the surface scaling relations (see the Appendix) and one-loop series for the bulk 
critical exponents obtained in Ref. [1, 2], 



1 v fi(a) — f 2 (a) 

4 2 

V = g/afa)™ + ... (55) 

The results of our calculation of [1/0] and [0/1] Pade-approximants of the series of the 
surface critical exponents at the adsorption threshold, and a group of critical exponents 
connected with the crossover exponent <3> are presented in Table 1 and Table 2, respectively. 

In the case a = d = 3, which corresponds to random uncorrelated point-like disorder 
(or short-range-correlated disorder) the obtained one-loop results for the surface critical 
exponents coincide with the results for the pure model (see [19, 24] Pade approximants 
[1/0] and [0/1]), as they should. For 2.3 < a < 3 the surface critical exponents, belonging to 
the new universality class associated with the LR fixed point of the RG equations, depend 
on a, similarly as the bulk exponents. This fact indicates that all the characteristics of the 
process of adsorption on a clean wall, described in the introduction, depend on the range of 
the correlations between the defects in the bulk. 



V. DISCUSSION OF THE RESULTS 



Let us first discuss the effect of the long-range correlated disorder on the distribution of 
monomers at and above the adsorption threshold (c > 0). From Eq.(5) and the fact that 
v{a) is a decreasing function [1, 2], it follows that the polymer with one end attached to 
the surface swells as in the bulk when a decreases. The behavior of the average end-to-end 
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distance at and above the threshold is independent of the surface critical exponents (see Eq. 
(5)), and the relevant exponent v(a) has been obtained at two loop order [1, 2]. 

Right at the threshold the average number of polymer ends and the average number 
of monomers depend on the distance from the surface. The distribution of monomers at 
different distances from the surface at the adsorption threshold, as well as the crossover 
behavior to the adsorbed state, can only be determined with the help of the surface critical 
exponents, as discussed in some detail in the introduction. The latter exponents have been 
calculated in this work up to one loop order and we can describe the effect of the long-range 
correlated disorder on the basis of our results. Our one-loop results show that for decreasing 
a the exponent 71(a) strongly increases, whereas $(a) decreases. We should point out that 
we cannot exclude the possibility that the dependence of the surface critical exponents on 
a is different beyond the one-loop approximation. 

Let us describe the effect of the long-range correlated disorder at the adsorption threshold 
and at the crossover to the adsorbed state assuming that the qualitative trends are properly 
captured by the one-loop results. For small distances (I <C z « £ R ) and for c = the 
partition function Zn(z) and the number of monomers in the layer at z are shown in Figs. 
1 and 2 respectively for several values of a. From these plots (see also Eqs. (7) and (8) and 
Tables I, II and Table in Ref.[l]) we can see that the number of free ends and the number of 
monomers in the near-surface region (/ <C z <C £r) both increase for deceasing a. Moreover, 
the dependence on z is the stronger, the larger the range of correlations between the defects. 
This result suggests that the larger is the range of correlations between the defects, the more 
efficient is the trapping of the chain between the attractive surface and the region occupied 
by the defects. For smaller a more steps are necessary for the chain to go around the 
defected region, which is larger for smaller a. From the fact that $ decreases for decreasing 
a it follows that the fraction of the monomers adsorbed at the surface also decreases (see 
Eq.(6) and Table II). This result is somewhat surprising, since it shows a trend just opposite 
to the one for the fraction of monomers contained in the layer at the distance z from the 
wall for I <C z <C £r. A possible explanation might be that once the polymer leaves the 
surface, the large disordered patches in the bulk make it difficult for the polymer to return 
back. As a result, the fraction of the monomers near the wall can be higher than right at 
the wall. Recall that also in the pure systems for c = there is a discrepancy between the 
behavior of the fraction of monomers right at the surface, Ni/N ~ iV* -1 , and the fraction of 
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monomers at the distance from the surface of the order of the effective monomer dimension 
/, M^(l)/N ~ jy*-i+7i-i_ Since 71 > 1 beyond MF, the concentration of monomers close 
to the surface is larger than just at the boundary, and Aq/M^(/) ~ jV 1 ^ 71 — > for iV — > 00 
and c = 0. Our one- loop results indicate that this effect is enhanced in the presence of 
long-range correlated disorder. 

Let us now consider the crossover to the adsorbed state. Assume that one end of the 
infinite polymer chain is attached to the surface, and the temperature is decreased below the 
threshold, i.e. — 1 < c < 0. The thickness of the adsorbed layer (Eq.(4)) starts to decrease 
from infinity when |c| increases from zero (see Fig. 3). From Eq.(2) and Table II we see that 
for a given c such that |c| < 1, the thickness of the near-surface, polymer layer increases 
when the range of correlations increases (a decreases). Moreover, the dependence on a of 
the exponent v(a)/&(a) is stronger than in the case of v{a). Hence, the dependence of the 
thickness of the adsorbed layer £ ~ | c | on a is stronger than the corresponding 
dependence of the mean end-to-end distance ~ N"( a ) [ n the bulk. It indicates that the 
effect of the long-range disorder on the adsorbed layer just below the threshold is stronger 
than the corresponding effect on the polymer coil in the bulk. The fraction of monomers, on 
the other hand, decreases for a fixed temperature for decreasing a, because for decreasing a, 
$(a) decreases (see Eq. (6) and Table II). 

In Fig. 4 the dependence of the fraction of monomers on the thickness of the adsorbed 
layer just below the threshold is shown for two values of a. Note that for a given value of 
the fraction of monomers at the surface the thickness of the adsorbed layer increases for 
decreasing a. Let us compare two systems, characterized by different ranges of disorder- 
correlations, ai < 02- Each system contains a chain with N monomers. Assume finally 
that the same number of monomers Aq is adsorbed at the surface in each system (of course 
temperatures in these systems are different, see (6) and Table II). The number of monomers 
inside the solution, N — Aq is the same in the two considered cases. Since ai < a 2 , larger 
correlated patches have to be avoided by the chain in the first system. By analogy with the 
polymer properties in the bulk [1] it is natural to expect that the same number of monomers 
contained in the solution, N — Aq, must effectively occupy larger space in the case where 
larger correlated patches have to be surrounded by the chain. Our results are thus consistent 
with the former results for the bulk systems [1]. 

We conclude that the long-range correlated disorder has a significant effect on the ad- 
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sorption of polymers on the surface at and near the adsorption threshold. When one end of 
the polymer is attached to the surface, the perpendicular part of the average distance to the 
other end increases for increasing range of correlations between the disorder. At the same 
time, the fraction of monomers at the surface decreases for fixed temperature both at and 
below the threshold. Moreover, just at the threshold the monomer concentration near the 
wall (i.e. for I <C z <C £r) increases. We should stress that the above conclusions follow 
from the results obtained in the one-loop approximation, and the possibility that the real 
dependence of the surface critical exponents on a is different cannot be excluded. Two-loop 
calculations and/or computer simulations, going beyond the scope of this work, might help 
to draw definite conclusions. 



Appendix 

The individual RG series expansions for the other critical exponents can be derived 
through standard surface scaling relations [14] with d = 3 

77 + ?7|| 





v*- 


-2 + 77,1), 


111 = 


= 1/(1- 


-^11). 


7i = 


1/(2 - 


■v±), 


Ai = 




-ml 




A 


d+2-r] 


5i = 


Ji = 


d - 2 + 77|| ' 




Ax 


d — 77|| 


6u = 




d - 2 + 7711 



(A2.1) 



Each of these critical exponents characterizes certain properties of the semi-infinite sys- 
tems with long-range quenched disorder, in the vicinity of the critical point. The values v, 
77, and A = u(d + 2 — 77) / 2 are the standard bulk exponents. 



Acknowledgments 

This work was supported by NATO Science Fellowships National Administration under 
Grant No. 14/B/02/PL. One of us (AC) acknowledges a partial support by the Polish 

19 



Ministry of Sciences through the Research Project No. 4T09A06622. 



[1] V.Blavats'ka, C.von Feber, and Yu.Holovatch, J.Mol.Liq. 92, 77 (2001). 
[2] V.Blavats'ka,C.von Feber, Yu.Holovatch, Phys.Rev.E, 64, 041102 (2001). 
[3] P.-G. de Gennes, Scaling Concepts in Polymer Physics (Cornell University) Press, Ithaca, NY, 
1979). 

[4] E.Eisenriegler, K. Kremer, and K. Binder, J. Chem. Phys. 77, 12 (1982). 
[5] E.Eisenriegler, J.Chem.Phys. 79, 1052 (1983) 

[6] E.Eisenriegler, Polymers near surfaces, World Scientific Publishing Co.Pte.Ltd., 1993. 

[7] R.Hegger and p.Grassberger, J.Phys.A 27, 4069 (1994). 

[8] Y.Singh, S.Kumar, and D.Giri, J.Phys.A 32, L407 (1999). 

[9] Y.Singh, D.Giri, and S.Kumar, J.Phys.A 34, L67 (2001). 
[10] R.Rajesh, D.Dhar, D.Giri, S.Kumar, and Y.Singh, Phys.Rev.E 65, 056124 (2002). 
[11] D.Zhao and T.Lookman, K.De'Bell, Phys.Rev.A 42, 4591 (1990). 
[12] H.W.Diehl and S.Dietrich, Z.Phys.B 42, 65 (1981). 
[13] H.W.Diehl and S.Dietrich, Z.Phys.B 50, 117 (1983). 

[14] H. W. Diehl, in Phase Transitions and Critical Phenomena, edited by C. Domb and J. L. 

Lebowitz (Academic Press, London, 1986), Vol. 10, pp. 75-267. 
[15] H. W. Diehl and A. Niisser, Z. Phys. B 79, 69 (1990). 

[16] J.des Cloizeaux and G.Jannink, Polymers in Solution (Clarendon Press, Oxford, 1990); 

L.Schafer L, Universal Properties of Polymer Solutions as Explained by the Renormalization 

Group (Springer, Berlin, 1999). 
[17] PG.de Gennes, J.Phys. (Paris) 37, 1445 (1976); M.Daoud and PG. de Gennes, ibid. 38, 85 

(1977). 

[18] M.N.Barber, A.S.Guttman, K.M.Middlemiss, G.M.Torrie, and S.G.Whittington, J.Phys.A 11, 
1833 (1978). 

[19] H. W. Diehl and M. Shpot, Nucl. Phys. B 528, 595 (1998). 
[20] Y.Kim, J.Phys.C 16, 1345 (1983). 

[21] A. B. Harris, Z. Phys. B: Condens.Matter 49, 347 (1983). 
[22] A. B. Harris, Journ. Phys. C 7, 1671 (1974). 

20 



[23] R.Guida and J.Zinn Justin J.Phys.A 31, 8104 (1998). 

[24] M. Shpot, Cond. Mat. Phys. N 10, 143 (1997). 

[25] Z.Usatenko, Chin-Kun Hu, Phys.Rev.E 65, 056102 (2001). 

[26] V.V.Prudnikov, P.V.Prudnikov, and A.A.Fedorenko, J.Phys.A 32, L399 (1999); 

V.V.Prudnikov, P.V.Prudnikov, and A.A.Fedorenko, ibid. 32,8587 (1999); V.V.Prudnikov, 

P.V.Prudnikov, and A.A.Fedorenko, Phys.Rev.B 62, 8777 (2000). 
[27] J.Zinn- Justin, Euclidean Field Theory and Critical Phenomena (Oxford Univ. Press, New 

York, 1989). 

[28] C.Bagnuls and C.Bervillier, Phys.Rev.B 24, 1226 (1981). 

[29] A.Ciach and U.Ritschel, Nucl.Phys.B 489, 653 (1997). 

[30] A. Weinrib and B. I. Halperin, Phys. Rev. B 27, 413 (1983). 

Figures 




FIG. 1: The partition function Zn{z) just at the adsorption thereshold c = and for N = 100, as 
a function of z/l for 1 <C z/l <C N u and for different values of a. Zn(z) is dimensionless and I is 
the microscopic length scale. 
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FIG. 2: The density of monomers M x (z) in the layer at the distance z from the surface to which 
one end of the chain is attached for l<z/I< N u just at the threshold c = and for N = 100 for 
different values of a. M x (z) is in arbitrary units and z/l is dimensionless. 




FIG. 3: The dependence of the thickness of the adsorbed layer £/l on |c|, for c < 0, where 
c (x (T — T a )/T a is the reduced temperature distance from the threshold, for different values of a. 
Both quantities are dimensionless. 
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FIG. 4: The dependence of the fraction of monomers at the surface N\/N on the thickness of the 
adsorbed layer £/l for c < (i.e. below the threshold) and for two values of a, a = 3.0 and a = 2.3. 
Both quantities are dimensionless. 



Tables 



TABLE I: Surface critical exponents of the long-flexible polymer at the special transition c = c a 
for d = 3 up to one-loop order calculated at the pure (the case a = 3 with (v* = 1.632, w* = 0)) 
and LR stable fixed point for different fixed values of the correlation parameter, 2 < a < 3. 



a v* w* 


V\\ 




V± 








7n 


7i 




[I/O] 


[o/i] 


[I/O] 


[o/i] 


[I/O] 


[o/i] 


[1/0] 


[1/0] 


3.0 1.63 0.00 


-0.204 


-0.169 


-0.102 


-0.114 


0.250 


0.250 


0.704 


1.255 


2.9 4.13 1.47 


-0.342 


-0.255 


-0.171 


-0.146 


0.247 


0.248 


0.837 


1.418 


2.8 4.73 1.68 


-0.402 


-0.287 


-0.200 


-0.167 


0.244 


0.244 


0.891 


1.480 


2.7 5.31 1.81 


-0.468 


-0.319 


-0.233 


-0.189 


0.241 


0.241 


0.951 


1.550 


2.6 5.89 1.87 


-0.542 


-0.351 


-0.270 


-0.212 


0.238 


0.238 


1.018 


1.630 


2.5 6.48 1.89 


-0.620 


-0.383 


-0.308 


-0.236 


0.235 


0.236 


1.090 


1.715 



2.4 7.10 1.87 -0.704 -0.413 -0.350 -0.259 0.233 0.233 1.169 1.810 
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2.3 7.76 1.84 -0.793 -0.442 -0.394 -0.283 0.230 0.231 1.253 1.911 



TABLE II: Surface critical exponents of the long-flexible polymer at the special transition c = c a 
for d = 3 involving the RG function r/ g , calculated at the pure (the case a = 3) and LR stable fixed 
point for different fixed values of the correlation parameter, 2 < a < 3. 



a v* 


w* 


Vc 




Ql 






a \\ 




$ 




(1-*) 


1* 






[1/0] 


[o/i] 


[I/O] 


[o/i] 


[ 


1/0] 


[o/i] 


[1/0] 


[o/i] 


[1/0] 


[o/i] 


3.0 1.63 


0.00 


-0.362 


-0.266 


0.217 


0.280 


-0. 


.362 


-0.266 


0.421 


0.427 


0.954 


0.956 


2.9 4.13 


1.47 


-0.607 


-0.378 


0.031 


0.181 


-0. 


.607 


-0.378 


0.363 


0.379 


0.942 


0.945 


2.8 4.73 


1.68 


-0.710 


-0.415 


-0.045 


0.147 


-0. 


.710 


-0.415 


0.335 


0.358 


0.950 


0.953 


2.7 5.31 


1.81 


-0.822 


-0.451 


-0.128 


0.114 


-0. 


,822 


-0.451 


0.306 


0.337 


0.955 


0.957 


2.6 5.89 


1.87 


-0.943 


-0.485 


-0.219 


0.082 


-0. 


,943 


-0.485 


0.276 


0.317 


0.954 


0.956 


2.5 6.48 


1.89 


-1.066 


-0.516 


-0.314 


0.051 


-1 


,066 


-0.516 


0.247 


0.298 


0.944 


0.947 


2.4 7.10 


1.87 


-1.192 


-0.544 


-0.413 


0.023 


-1. 


.192 


-0.544 


0.222 


0.282 


0.922 


0.928 


2.3 7.76 


1.84 


-1.309 


-0.567 


-0.511 


-0.003 


-1 


.309 


-0.567 


0.203 


0.271 


0.881 


0.894 
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